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Abstract: We calculate the thermodynamic functions of pure-glue QCD to three- loop 
rS I order using the hard-thermal-loop perturbation theory (HTLpt) reorganization of finite 

C^ I temperature quantum field theory. We show that at three-loop order hard-thermal-loop 

perturbation theory is compatible with lattice results for the pressure, energy density, 
and entropy down to temperatures T ~ 3 Tc. Our results suggest that HTLpt provides 
a systematic framework that can be used to calculate static and dynamic quantities for 
temperatures relevant at LHC. 
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The goal of ultrarelativistic heavy-ion collision experiments is to generate energy den- 
sities and temperatures high enough to create a quark-gluon plasma. One of the chief 
theoretical questions which has emerged in this area is whether it is more appropriate to 
describe the state of matter generated during these collisions using weakly-coupled quan- 
tum field theory or a strong-coupling approach based on the AdS/CFT correspondence. 
Early data from the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Labs 
indicated that the state of matter created there behaved more like a fluid than a plasma 
and that this "quark-gluon fluid" is strongly coupled [||]. 

In the intervening years, however, work on the perturbative side has shown that observ- 
ables like jet quenching [^ and elliptic flow [?] can also be described using a perturbative 
formalism. Since in phenomenological applications predictions are complicated by the mod- 
eling required to describe, for example, initial-state effects, the space-time evolution of the 
plasma, and hadronization of the plasma, there are significant theoretical uncertainties re- 
maining. Therefore, one is hard put to conclude whether the plasma is strongly or weakly 
coupled based solely on RHIC data. To have a cleaner testing ground one can compare 
theoretical calculations with results from lattice quantum chromodynamics (QCD). 

Looking forward to the upcoming heavy-ion experiments scheduled to take place at the 
Large Hadron Collider (LHC) at the European Organization for Nuclear Research (CERN) 



o 
o 
o 



it is important to know if, at the higher temperatures generated, one expects a strongly- 
coupled (liquid) or weakly-coupled (plasma) description to be more appropriate. At RHIC, 
initial temperatures on the order of one to two times the QCD critical temperature, Tc ~ 
190 MeV, were obtained. At LHC, initial temperatures on the order of 4— 5 Tc are expected. 
The key question is, will the generated matter behave more like a plasma of quasiparticles 
at these higher temperatures. 

The calculation of thermodynamic functions using weakly-coupled quantum field the- 
ory has a long history |Q, |5|, ^, |^, ^, ^ ^, |T^, [l^, |l^, 14]. The QCD free energy is known up 
to order g^ log{g); however, the resulting weak-coupling approximations do not converge at 
phenomenologically relevant couplings. For example, simply comparing the magnitude of 
low-order contributions to the QCD free energy with three quark fiavors {Nf = 3), one finds 
that the ^f^ contribution is smaller than the g^ contribution only for g -^ 0.9 (a^ ~ 0.07) 
which corresponds to a temperature of T ~ 10^ GeV ~ 5 x 10^ Tc- 

In Fig 1^, we show the weak-coupling expansion for the pressure of pure-glue QCD 

5/2 

normalized to that of an ideal gas through order as ■ The various approximations oscillate 
wildy and show no signs of convergence in the temperature range shown. The bands are 
obtained by varying the renormalization scale ^u by a factor of two around the value /i = 2ttT 
and we use three-loop running of a^ This oscillating behavior is generic for hot field theories 
and not specific to QCD. 
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Figure 1: Weak-coupling expansion for the scaled pressure of pure-glue QCD. Shaded bands show 
the result of varying the renormalization scale /x by a factor of two around fj, = 27rT. 



The poor convergence of finite-temperature perturbative expansions of thermodynamic 
functions stems from the fact that at high temperature the classical solution is not described 
by massless gluons. Instead one must include plasma effects such as the screening of electric 
fields and Landau damping of excitations via a self-consistent hard-thermal-loop (HTL) 



resummation |15]. There are several ways of systematically reorganizing the perturbative 
expansion |]l^] . Here we will present the details of a new NNLO calculation which uses the 
hard-thermal- loop perturbation theory method |17, 18, |l^, 22 1 and compare with previously 



obtained LO and NLO results. 

The basic idea of the technique is to add and subtract an effective mass term from the 
bare Lagrangian and to associate the added piece with the free part of the Lagrangian and 
the subtracted piece with the interactions [^, 24 1. However, in gauge theories, one cannot 



simply add and subtract a local mass term since this would violate gauge invariance |25 
|2q , P7| ]. Instead, one adds and subtracts an HTL improvement term which modifies the 
propagators and vertices self-consistently so that the reorganization is manifestly gauge 



invariant [25|. 



1. HTL perturbation theory 

The Lagrangian density for pure-glue QCD in Minkowski space is 



£gf = --Tr 



{d^^Af^Y . (1.2) 



The perturbative expansion in powers of g generates ultraviolet divergences. The 
renormalizability of perturbative QCD guarantees that all divergences in physical quantities 
can be removed by renormalization of the coupling constant a^ = g^/47r. There is no 
need for wavefunction renormalization, because physical quantities are independent of the 
normalization of the field. There is also no need for renormalization of the gauge parameter, 
because physical quantities are independent of the gauge parameter. 

Hard-thermal-loop perturbation theory (HTLpt) is a reorganization of the perturbation 
series for thermal QCD. The Lagrangian density is written as 

^ = (>Cqcd + >Chtl) ^ +A/:htl. (1.3) 
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In this paper we discuss the calculation of thermodynamic functions of a gas of glu- 
ons at phenomenologically relevant temperatures using hard-thermal-loop perturbation 
theory. We present results at leading order (LO), next-to-leading order (NLO), and next- ' \J 

to-next-to-leading order (NNLO) and compare with available lattice data [^ ^ for the (^ ^ 

thermodynamic functions of SU(3) Yang- Mills theory. The calculation is based on a reorga- 
nization of the theory around hard-thermal-loop (HTL) quasiparticles. Our results indicate 
that the lattice data at temperatures T ~ 2 — 3 Tc are consistent with the quasiparticle 
picture. This is a non-trivial result since, in this temperature regime, the QCD coupling L^-' 

constant is neither infinitesimally weak nor infinitely strong with 5 ~ 2, or equivalently C !) 

as = (7^/(47r) ~ 0.3. Therefore, we have a crucial test of the quasiparticle picture in the 
intermediate coupling regime. 
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£qcd = -^TV [G^.G^1 + /Igf + /Igh + A/:qcd , (1-1) O 

where the field strength is G^^'^ = d^^A^ — d'^A'^ — ^g[A^^,A'^]. The ghost term /3gh depends 
on the gauge- fixing term £gf . In this paper we choose the class of covariant gauges where 
the gauge-fixing term is 



The HTL improvement term is 



CuTL = -kl - 5)mlTr ( G^„ (7^^ ) G% 1 , (1.4) 




where the covariant derivative is D^ = d^ — igA^ and y'^ = (1, y) is a hght-like four- vector, 
and {■ ■ ■)y represents the average over the directions of y. The term (1^) has the form 



of the effective Lagrangian that would be induced by a rotationally invariant ensemble of 
charged sources with infinitely high momentum. The parameter mo can be identified with 
the Debye screening mass. HTLpt is defined by treating 5 as a formal expansion parameter. 
The HTL perturbation expansion generates ultraviolet divergences. In perturbative 
QCD, renormalizability constrains the ultraviolet divergences to have a form that can be 
cancelled by the counterterm Lagrangian A/Iqcd- We will demonstrate that renormalized 
perturbation theory can be implemented by including a counterterm Lagrangian A>Chtl 



M 

among the interaction terms in (L3). There is no proof that the HTL perturbation expan- ' \J 

sion is renormalizable, so the general structure of the ultraviolet divergences is not known; ^~^ 

however, it was shown in previous papers |17, |l8| that it was possible to renormalize the ^ — >. 
NLO order HTLpt prediction for the free energy of QCD using only a vacuum counterterm, 

a Debye mass counterterm, and a fermion mass counterterm. In this paper we will show """^ 

that this is also possible at NNLO. In particular, the only new counterterm we need to L\J 

introduce is for the coupling constant as , which coincides with its perturbative value giving C !) 
rise to the standard one-loop running. 

We find that the counterterms necessary to renormalize HTLpt at NNLO are 

(JAas = -^^ald"^ + 0{5^al) , (1.5) ' 

Izvre 

Ami, = (-^«^^ + ^('^'«')) (1 - ^Wd , (1-6) Q 



A£-o = ^y^ + 0{6as) ] (1 - 6rm% . (1.7) 



We note that the counterterm in Eq. ( |1.5[ ) coincides with the perturbative one-loop running. 

Physical observables are calculated in HTLpt by expanding them in powers of 5, trun- 
cating at some specified order, and then setting 5 = 1. This defines a reorganization of 
the perturbation series in which the effects of mj^ term in (1.4) are included to all orders 
but then systematically subtracted out at higher orders in perturbation theory by the Srnjj 
terms in (|l.4D . If we set (5 = 1, the Lagrangian (|l.3D reduces to the QCD Lagrangian (|L1|). 

If the expansion in 6 could be calculated to all orders the final result would not depend 
on m£) when we set 6 = 1. However, any truncation of the expansion in 5 produces results 
that depend on rriD ■ Some prescription is required to determine ttid as a function of T and 
Og. We will discuss several prescriptions in Sec. VI. 

2. Diagrams for the thermodynamic potential 

In this section, we list the expressions for the diagrams that contribute to the thermody- 
namic potential through order 5'^ in HTL perturbation theory. The diagrams are shown 



o 



in Figs. 1^, and y. A key to the diagrams is given in Fig. ^. The expressions here wih be 
given in Euclidean space; however, in Appendix ^ we present the HTLpt Feynman rules 
in Minkowski space. 
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Figure 2: Diagrams contributing through NLO in HTLpt. The spiral lines are gluon propagators 
and the dotted lines are ghost propagators. A circle with a 11 indicates a gluon self-energy insertion 
All propagators and vertices shown are HTL-resummed propagators and vertices 
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Figure 3: Diagrams contributing to NNLO in HTLpt which contribute through order g^. The 
spiral lines are gluon propagators and the dotted lines are ghost propagators. A circle with a H 
indicates a gluon self-energy insertion The propagators are HTL-resummed propagators and the 
black dots indicate HTL-resummed vertices. The lettered vertices indicate that only the HTL 
correction is included. The yellow box denotes the insertion of the one-loop self-energy defined in 
Fig. I 
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Figure 4: Key to the diagrams in Figs, g and |[ 
The thermodynamic potential at leading order in HTL perturbation theory for QCD 
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^Lo = {Nl - 1) J-ia+ib + Ao£-o . (2.1) 

Here, J^ia+ib is the contribution from the gluon and ghost diagrams shown on the first hne 
of Fig. I 

-^ia+ib = -^^ {{d - 1) log [-At(P)] + log Al(P)} . (2.2) 

The transverse and longitudinal HTL propagators Aj'(P) and Al{P) are given in ( |A.4g| ) 
and ([A. 50). The leading-order vacuum counterterm AqSq was determined in Ref. |J17|: 
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The thermodynamic potential at NLO in HTL perturbation theory can be written as [-i— | 

+Aifo + Aiml-^n^o , (2.4) O 

where Ai£'o and Aim'j-, are the terms of order 5 in the vacuum energy density and mass [S^ 

o 
I— ^ 
o 



counterterms: 



A 2 ll^c 2 / X 

Aim^ = --- — agmjy . (2.6) 

izTre 



iV. 



^^ 



^2b = ^-^g't r^^'^'^^P, -P, Q, -Q)A^-(P) 

xA^"(Q) , (2.8) 

where R = Q — P. For the corresponding diagrams, see the second line of Fig. |2|. 
The contribution from the ghost diagram is 

^2c = ^9'f^ ^j^Q^R-'A^^'^iP). (2.9) 

The contribution from the HTL gluon counterterm diagram with a single gluon self- 
energy insertion is 

•^2d = l'^ n^^(P)A^'^(P) . (2.10) 



O 

The contributions from the two-loop diagrams with the three-gluon and four-gluon vertices (^ ^ 

are 

J-^ TPQ 

xA^''{Q)AP^R) , (2.7) 



o 



The thermodynamic potential at NNLO in HTL perturbation theory can be written 
as 

^NNLO = f^NLO + (^c ~ 1) [-^Sa + -^3b + -^3c + -^3d + -^3c + -^3f + •^3g + -^3h + -^3i + -^3j + -^3k 

d d 



1 / 52 



,a^o) {^imlf + {Nl - l)^^^±?l^Aia. . (2.11) 



2 \{dm]jY^ 

where Aiq^, A2m|), and /S.2^0 are the terms of order 5"^ in the couphng constant, mass, _ 

and vacuum energy density counterterms: _^ j 

^2ml = ^^asml , (2.13) 'Xl 

A^2 _ 1 , O 

^=^° = i4a"* ■ <^ "' o 

The contributions from the three-loop diagrams are given by ^,,_^ 

NJ 

-^3a = ^g% T'^^HP, Q, -P - Q)A"^(P)A^^(Q)A^-(P + Q)T^'-\-Q, -R, Q + R)A^-{R) ^ 
^4 TpQR 

xA^^iQ + R)T''^P{P + Q,-Q-R,R- P)A^'^(i? - P)T^'t"'{-P, P-R,R) , (2.15) 



O 



xA''^(P + Q)A^^(Q), (2.16) O 

-^3d = T^X r°^'^^(P,Q,i?,S)rT^''"^(P,Q,i?,5)A"T(P)A^^(Q)A^'^(i?)A'^^(5) , (2.18) 

A^2_f 
J-ge = — f X r"^'^'^(P, g, P, 5)A"^(P)A^'^(Q)A^'5(P)A<^'^(5)A^^(P + Q) 

xr/5'^^(_P, _Q, p + Q)r^'^'^(-P - Q, -P, -S) , (2.19) 

7:3^ = 4 n^^(p)A^°(p)n"''(p)A^'^(p) , (2.20) 

-^3g--^5^^^^ p4(p + g)2(P + p)2 ^ (^)^ (^)' (2.21) 

where S = —{P + Q + R) and n^'^(P) is the one- loop gluon self-energy defined by the 
yellow box in Fig. ^ 

n^-(P) = \N,g^i r^^'^<-(^{P,-P,Q,-Q)A"f'{Q) + ^N,g^4 T'^'^^ {P, Q , - P - Q) A""' (Q) 
xr'^^'^(P, Q, -P - Q)AH{R) + N,g'^ q2[pIq]1 ' (^.22) 



The contributions from the two-loop diagrams with a single self-energy insertion are 

4 TPQ 

J-g; = -^5'rf r"'^'^'^(P,-P,Q,-Q)A"^(P)nT^(P)A^^(P)A^''(Q) , (2.24) 

4 TpQ 

^'=„2^c/ P"(P + Q) AaM//^^TTMl'//^^ Ai'/?/ 






-^3j = -^rj^ p2(p + QJ2 ^"^(Q)n^'(Q)^'^(Q) ' (2-25) 

where R = Q — P. 

The two-loop diagrams with a subtracted vertex is C__| 

xA^'^(Q)A^^(P) , (2.26) 'Xl 

7-31 = ^g^mli r^-'^-(P, -P, Q, -Q)A'^''(P) O 

« TpQ Q 

xA^'^(Q) , (2.27) ^_^ 

where R = Q-P. ^ 

The contribution from the HTL gluon counterterm diagram with two gluon self-energy C) 

insertions is 

•^3m = -\i n^'^(P)A^°(P)n"'^(P)A'^'^(P) . (2.28) 

3. Expansion in rriij/T 



was clear that evaluating these scalar sum-integrals exactly was intractable and the sum- 
integrals were calculated approximately by expanding them in powers of m£)/T. We will 
follow the same strategy in this paper and carry out the expansion to high enough order 
to include all terms through order g^ if mD is taken to be of order g. 

The free energy can be divided into contributions from hard and soft momenta. In 
the one- loop diagrams, the contributions are either hard [h) or soft (s), while at the two- 
loop level, there are hard-hard {hh), hard-soft (/is), and soft-soft {ss) contributions. At 
three loops there are hard-hard-hard (hhh), hard-hard-soft (hhs), hard-soft-soft (hss), and 
soft-soft-soft (sss) contributions. 

3.1 Leading order 

3.1.1 Hard contribution 

For hard momenta, the self-energies are suppressed by m£)/T relative to the propagators, 
so we can expand in powers of IiT{P) and n£,(P). 

For the one-loop graphs (la) and (lb), we need to expand to second order in m?j^: 
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o 
o 

In the papers [17, 18, 19, 20, 21|, the free energy was reduced to scalar sum-integrals. It C !) 
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\{d-l)^^\og{P^ 
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p4 



,2p2 
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r)2p2 



p- 



l+|2 + 2|(zi)u 



, 2e 



\ATrTJ ^ 



2e 



--7 + 27. + ^) (^)' mi, 



(3.1) 



3.1.2 Soft contribution 



The soft contribution in the diagrams (la + 16) arises from the Pq = term in the sum- 
integral. At soft momentum P = (0, p), the HTL self-energy functions reduce to HxiP) = 



and Ul{P) 



m 



D- 



The transverse term vanishes in dimensional regularization because 



there is no momentum scale in the integral over p. Thus the soft contributions come from 
the longitudinal term only and read 



-F 



is) 



la+lb 



h: j\og{p' + ml) 



m-^jjT 



Uir \2mJ 



(J^\ 



2e 
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(3.2) 



We have kept the order-e in the m'j-j and mj-) terms, respectively in Eqs. ( ^.1| ) and 
since they contribute in the counterterms at next-to-leading order. 

3.2 Next-to-leading order 

3.2.1 Hard contribution 

The one-loop graph with a gluon self-energy insertion (2d) has an explicit factor of mj, 
and so we need only to expand the sum-integal to first order in m'j-,: 



o 
o 

o 



O 

o 
o 
o 



J"' 



(h). 



\f 1 



2d 



-2'^%P^^W 



1 



4 

-mj, 



24 [ V C(-l) 

3.2.2 Soft contribution 
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(4^^""™"^^ + 
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647r2 



p2p2 

--7 + 2-iE + 



1 - ^ ?^U- 



/^ 



KAt^TJ 



T^S') 



The soft contribution from (2d) arises from the Pq = term in the sum-integral. Only the 
longitudinal part n^, (P) of the self-energy contributes and reads 
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Stt \2mDj ^ ^ 



(3.4) 
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3.2.3 {hh) contributions 

For hard momenta, the self-energies are suppressed by ruD/T relative to the propagators, 
so we can expand in powers of 11^ and 11^. The two- loop contribution was calculated in 
Ref. ilHl and reads 
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2a+2b+2c 
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p2Q2^4 



+ 4 



p2g2^2^2 
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1 1 
P^lPQ^ 
1 1 



Tq 



Tr +4d „,L ^ Tfl + 2d-^-^TR 



p2Q2^2'^ p2Q2j.4'^ p2Q2j.4 

Using the expressions for the sum- integrals listed in Appendix B, we obtain 



P2 q2Q2 

(3.5) 
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37r V47rr/ 
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(3.6) 



3.2.4 (/is) contributions 

In the (hs) region, the momentum P is soft. The momenta Q and R are always hard. 
The function that multiplies the soft propagator At(0, p), Ai(0,p) or Ax(0, p) can be 
expanded in powers of the soft momentum p. In the case of At (0, p) , the resulting integrals 
over p have no scale and they vanish in dimensional regularization. The integration measure 
J scales like m^, the soft propagators Ai(0, p) and Ax(0, p) scale like l/mj-,, and every 
power of p in the numerator scales like m£) . The two-loop contribution was calculated in 



Ref. 1 18] and reads 
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~T9 ^ / 2 , 2 



-(^-l)^ + 2(d-l)|;' 



Ncg^mlT 
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~d Q8 



(3.7) 



In order to facilitate the calculations, it proves useful to isolate the terms that are 
specific to HTL perturbation theory. After integrating by parts and using the results from 
Zhai and Kastening [|ll|, we can write 
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2a+2b+2c 



^^'^^'-'^'L^A 



Q2 + 12 ^"^ 



-%d-5)g^Tml[^^i^^. 
2 JpP^ + mj^TQ Q^ 



(3.8) 



Using the expressions for the integrals and sum-integrals in Appendices B and C, we obtain 
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3.2.5 (ss) contribution 



The (ss) contribution to tlie free energy is given by a two-loop calculation in electrostatic 
QCD (EQCD) in three dimensions. This calculation was carried out in Ref. Q by Braaten 
and Nieto. Alternatively, one can isolate the (ss) contributions from the two-loop diagrams 



which were calculated by Arnold and Zhai in Ref. H. In Ref. |18], this contribution was 
calculated and agrees with earlier results. One finds 
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We have kept the order e in terms Eqs. (3^), (|3.4|), (^^), and (|3^ 
in the counterterms at NNLO. 



(3.10) 
since they contribute 



3.3 Next-to-next-to-leading order 

3.3.1 Hard contribution 

The one-loop graph with two gluon self-energy insertions (3m) is proportional to Tiij^ and 
so must be expanded to zeroth order in mj^ 



Tl 



(h) 



1 



3m 



4(d- 
1 



1 



4 



1 1 

9 

p4 p2p2 



1 „ 2tt^ 

- - 7 + 27ij + — - 
e 3 



2'^^'7> + 2^-7'p + rf^ (Tpf 



(3.11) 



1287r2 

3.3.2 Soft contribution 

The soft contribution from (3m) arises from the Pq = term in the sum-integral. Only the 
longitudinal part liLiP) of the self-energy contributes and reads 
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(3.12) 



3.3.3 {hh) contributions 

We also need the (hh) contribution from the diagrams (3h) — (31). We calculate their 
contributions by expanding the two- loop diagrams (2a) — (2c) to first order in m^. This 
yields 



Tl 



(hh) 



N, 



3h-31 



c 2 2 

-9 m^ 



PQ 



-2(<i-l);^^+2(d-2)Jj^ + (d + 2; 



1 



Q2i?2^2 



-2d 



P-Q 

p2Q2^4 



4d 



p2Q2^4 



+ 4 



p2Q2j,2Jl2 



2{d- 1) 



1 1 

P2^2Q2 



Tq 



7 
96 



[d+l] 
"1 



P'^Q-r 






Tr + M 



+ 4.621 



Nm, 



37r V47rr/ 



2^jR + '^d 

4e 



P2g-r 



p2g2^4 



Tr 



—-= runT . 



D^ 



(3.13) 



o 
o 

o 



o 

o 
o 
o 



12 



3.3.4 (/is) contributions 

We also need the (hs) contribution from the diagrams (3h) — (31). Again we calculate their 
contributions by expanding the two- loop diagrams (2a) — (2c) to first order in m'j-,. This 
yields 



r. 
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2 Jp [p^ + mj^y 
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(3.14) 



This yields 
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3.3.5 (ss) contribution 

The (ss) contribution from the two-loop diagrams with a single self-energy insertion or 
vertex correction can be obtained by expanding the two-loop result in powers of rnj^. This 
yields 

4 2(p2 _^ 4^2^) 
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We have verified this by explicitly calculating the relevant diagrams. 
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3.3.6 (hhh) contribution 

If all the three loop momenta are hard, we can obtain the m£)/T expansion simply by 
expanding in powers of m'j^. To obtain the expansion through order g^, we can use bare 
propagators and vertices. The contributions from the three-loop diagrams were first cal- 
culated by Arnold and Zhai in Ref. [g], and later by Braaten and Nieto Q. One finds 
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3.3.7 {hhs) contributions 

All the diagrams except (3f) are infrared finite in the limit rriD — )> 0. This implies that 
the g^ contribution is given by using a dressed longitudinal propagator and bare vertices. 
The ring diagram (3f) is infrared divergent in that limit. The contribution through g^ is 
obtained by expanding in powers of self-energies and vertices and one obtains 
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3.3.8 {hss) contribution 

For all the diagrams that are infrared safe, the (hss) contribution is of order g'^m'jj, i.e. g' 
and can be ignored. The infrared divergent diagrams contribute as follows 
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3.3.9 {sss) contribution 

The (sss) contribution to the free energy is given by a three-loop calculation of the free en- 
ergy of Electrostatic QCD in three dimensions. This calculation was performed in Ref. Q. 
Alternatively, one can isolate the {sss) contributions from the diagrams listed in Ref. Q. 
The result is 
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The expression for the integrals are given in Appendix C. Adding Eqs. ( |C.ll )- ( |C.23| ), 
the final result is 
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Note that all the poles in e cancel. 
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4. The thermodynamic potential 

In this section we present the final renormalized thermodynamic potential explicitly through 
order 5^, aka NNLO. The final NNLO expression is completely analytic; however, there are 
some numerically determined constants which remain in the final expressions at NLO. 

4.1 Leading order 

The leading order thermodynamic potential is given by the contribution from the diagrams 
(la) and (lb). 
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(4.1) 

where Jideai is the free energy of a gas of A''^ — 1 massless spin-one bosons and ttid and jl 
are dimensionless variables: 
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The complete expression for the leading order thermodynamic potential is given by adding 

(4.5) 



the leading vacuum energy counterterm (p.3|) to Eq. (4.1): 
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4.2 Next-to-leading order 

The renormalization contributions at first order in 6 are 
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Using the results listed in Eqs. ( |2.5D and (|2.6D the complete contribution from the coun- 
terterm at first order in 6 is 
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The NLO thermodynamic potential reads 
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4.3 Next-to-next-to-leading order 

The renormalization contributions at second order in 6 are 
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Using the results listed above, we obtain 
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Adding the NNLO counterterms ( [4.10 ) to the contributions from the various NNLO 
diagrams we obtain the renormalized NNLO thermodynamic potential. We note that at 
NNLO all numerically determined coefficients of order e'^ drop out and we are left with a 
final result which is completely analytic. The resulting NNLO thermodynamic potential is 
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Note that if we use the weak-coupling value for the Debye mass ml, = AirNcasT'^/'i, 
the NNLO HTLpt result ( [4.11 ) is guaranteed to reduce to the weak-coupling result through 
order g^ and we have checked that this is the case. 

5. Thermodynamic functions 
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5.1 Mass prescriptions 

The mass parameter ttid in HTLpt is completely arbitrary. To complete a calculation, it 
is necessary to specify m/5 as function of g and T. In this section we will discuss several 
prescriptions for the mass parameter. 



5.1.1 Variational Debye mass 

The variational mass is given by the solution to the equation 
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At leading order in HTLpt, the only solution is the trivial solution, i.e. rriD = 0. In 
that case, it is natural to chose the weak-coupling result for the Debye mass. This was 
doneinRef. [|l7|. 

At NLO, the resulting gap equation has a nontrivial solution, which is real for all values 
of the coupling. At NNLO, the solution to the gap equation ( |5.2| ) is plagued by imaginary 
parts for all values of the coupling. The problem with complex solutions seems to be generic 
since it has also been observed in screened perturbation theory |24] and QED [19, pH, 21 



In those cases, however, it was complex only for small values of the coupling. 

5.1.2 Perturbative Debye mass 

At leading order in the coupling constant g, the Debye mass is given by the static longitu- 
dinal gluon self-energy at zero three-momentum, mj-, = 112,(0,0), i.e. 
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^N^asT^ . (5.3) ^ 



3 

The next-to-leading order correction to the Debye mass is determined by resummation of 
one- loop diagrams with dressed vertices. Furthermore, since it suffices to take into account 
static modes in the loops, the HTL-corrections to the vertices also vanish. The result, 
however, turns out to be logarithmically infrared divergent, which reflects the sensitivity 
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to the nonperturbative magnetic mass scale. The result was first obtained by Rebhan |43| 

and reads C ^ 

2m£) 1 



o 

where mmag is the nonperturbative magnetic mass. We will not use this mass prescription C ^ 

since it involves the magnetic mass which would require input from e.g. lattice simulations. ^~^ 

5.1.3 BN mass parameter m'^ 

In the previous subsection, we saw that the Debye mass is sensitive to the nonperturbative 
magnetic mass which is of order g^T. In QED, the situation is much better. The Debye 
mass can be calculated order by order in e using resummed perturbation theory. The 
Debye mass then receives contributions from the scale T and eT. Effective field theory 
methods and dimensional reduction can be conveniently used to calculate separately the 
contributions to mu from the two scales in the problem. The contributions to mu and 
other physical quantities from the scale T can be calculated using bare propagators and 
vertices. The contributions from the soft scale can be calculated using an effective three- 
dimensional field theory called electrostatic QED. The parameters of this effective theory 
are obtained by a matching procedure and encode the physics from the scale T. The 
effective field theory contains a massive field Aq that up to normalization can be identified 
with the zeroth component of the gauge field in QED. The mass parameter niE of ^o gives 
the contribution to the Debye mass from the hard scale T and can be written as a power 
series in e^. 
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For nonabelian gauge theories, the corresponding effective three-dimensional theory is 
cahed electrostatic QCD. The mass parameter niE for the field Aq (which lives in the adjoint 
representation) can also be calculated as a power series in g^. It has been determined to 
order g'^ by Braaten and Nieto B. For pure-glue QCD, it reads 
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We will use the mass parameter niE as another prescription for the Debye mass and denote 
it by the Braaten Nieto (BN) mass prescription. 

5.2 Pressure 

In this subsection, we present our results for the pressure using the variational mass pre- 
scription and the BN mass prescription. 

5.2.1 Variational mass 
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Figure 5: Comparison of LO, NLO, and NNLO predictions for the scaled real part of the pressure 
using the variational mass prescription. Shaded bands show the result of varying the renormalization 
scale /i by a factor of two around /i = I'kT . 

In Fig. ^, we compare the LO, NLO, and NNLO predictions for the real part of the 
pressure normalized to that of an ideal gas. Shaded bands show the result of varying the 
renormalization scale /i by a factor of two around [i = 2txT . We use three-loop running of 

In Fig. g, we show the NNLO result for the imaginary part of the pressure normalized 
by the ideal gas pressure. We use three- loop running of a^. The imaginary part decreases 
with increasing temperature and is rather small beyond 3-4 T^,- 
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Figure 6: The NNLO result for the scaled iraaginary part of the pressure with three- loop running 
and variational mass prescription. The two curves arise from the two complex conjugate solutions 
to the gap equations. Shaded bands show the result of varying the renormalization scale /i by a 
factor of two around /z = 2'kT. 
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Figure 7: Comparison of LO, NLO, and NNLO predictions for the scaled pressure using the BN 
mass prescription and one-loop running of as . The points are lattice data for pure-glue with Nc — 3 
from Boyd et al. [Q. Shaded bands show the result of varying the renormalization scale /i by a 
factor of two around ^ = 2'kT . 

5.2.2 BN mass 

In Fig. 0, we show the HTLpt predictions for the pressure normalized to that of an ideal gas 
using the BN mass prescription. The bands are obtained by varying the renormalization 
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scale by a factor of two around ^ = 2itT. We use one-loop running of a^. In Fig. P, we 
again plot the normalized pressure, but now with three-loop running of a^. The agreement 
between the lattice data from Boyd et al. [^] is very good down to temperatures of around 
3 Tc- Comparing Figs. |^-^ we see that using the three-loop running, the band becomes 
wider. However, the difference is significant only for low T, where the HTLpt results 
disagrees with the lattice anyway. For T > 3Tc, the prescription for the running makes 
very little difference. 
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Figure 8: Comparison of LO, NLO, and NNLO predictions for the scaled pressure using the BN 
mass prescription and three-loop running of as ■ The points are lattice data pure-glue with N^ = 3 
from Boyd et al. pSi . Shaded bands show the result of varying the renormalization scale /j, hj a. 
factor of two around /i — 2ttT. 

Until recently, lattice data for thermodynamic variables only existed for temperatures 
up to approximately 5 Tc. In the paper by Enrodi et al |2^, the authors calculate the 
pressure on the lattice for pure-glue QCD at very large temperatures. In Fig. ^, we show 
the results of Enrodi et al as well as Boyd et al, together with the HTLpt NLO and NNLO 
predictions for the pressure. The two points from Ref. p^ have large error bars, but the 
data points are consistent with the HTL predictions. 

We note that our prediction for the normalized free energy using either the leading- 
order, BN, or variational mass prescriptions is independent of Nc if one holds OgNc fixed 
('t Hooft limit). However, this need not be the case for an arbitrary mass prescription. 
The A'^c-independence of the normalized free energy of the free energy is in agreement 
with recent lattice measurements that show a very small dependence of the free energy on 
Nc ^,^. 

Due to the imaginary parts we do not show results coming from the variational pre- 
scription in the remainder of the results section. We will return to a discussion of the 
dependence on mass prescriptions in the conclusions. 
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Figure 9: Comparison of NLO, and NNLO predictions for the scaled pressure with SU(3) pure-glue 
lattice data from Boyd et al. pSJ and Endrodi et al. |M. Shaded bands show the result of varying 
the renormalization scale /i by a factor of two around ^ — 2'itT. 

5.3 Energy density 

The energy density £ is defined by 



£ = T -T—- 
dT 



(5.6) 



In Fig. 10, we show the LO, NLO, and NNLO predictions for energy density normalized 
to that of an ideal gas. We use three-loop running and the BN mass. The bands show 
the result of varying the renormalization scale /x by a factor of two around /i = lnT . Our 
NNLO predictions are in excellent agreement with the lattice data down to T ~ 2Tc. 

5.4 Entropy 

The entropy density is defined by 



S = -—- 



d7_ 
dT ' 



(5.7) 



where all other parameters that J- depends on, are kept fixed. In Fig. 11, we show the 



entropy density normalized to that of an ideal gas. We use three-loop running and the BN 
mass. The points are lattice data from Boyd et al. [g^]. Our NNLO predictions are in 
excellent agreement with the lattice data down to T ~ 2Tc. 

5.5 Trace anomaly 

In pure-glue QCD or in QCD with massless quarks, there is no mass scale in the Lagrangian 
and the theory is scale invariant. At the classical level, this implies that the trace of the 
energy-momentum tensor vanishes. At the quantum level, scale invariance is broken by 
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Figure 10: Comparison of LO, NLO, and NNLO predictions for the scaled energy density with 
SU(3) pure-glue lattice data from Boyd et al. 28 . We use three-loop running and the BN mass. 



Shaded bands show the result of varying the renormalization scale /i by a factor of two around 
M = 27rT. 
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Figure 11: Comparison of LO, NLO, and NNLO predictions for the scaled entropy with SU(3) 
pure-glue lattice data from Boyd et al. pa. We use three-loop running and the BN mass. Shaded 
bands show the result of varying the renormalization scale /i by a factor of two around /i = 2ttT. 



renormalization effects. It is convenient to introduce the scale anomaly density £ — 3V, 
which is proportional to the trace of the energy-momentum tensor. The trace anomaly can 
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be written as 



S-3V = -T^— ( — 



(5.8) 



In Fig. |1^, we show the HTLpt predictions for the trace anomaly divided by £^idcai using 
the BN mass prescription and three-loop running of a^. The points are lattice data from 
Boyd et al. [^] For temperatures below approximately 2Tc, there is a large discrepancy 
between the HTLpt predictions and lattice data. At LO and NLO, the curves are even 
bending downwards. 

At temperatures close to the phase transition it has been suggested that the discrepancy 
between HTLpt resummed predictions for thermodynamics functions and, in particular, 
the trace anomaly is due to influence of a power corrections ||3^ , 33, 34, |3^, ^ which are 
related to confinement. Phenomenological fits of lattice data which include such power 
corrections show that the agreement with lattice data is improved p^, p8|. Alternatively, 
others have constructed AdS/CFT inspired models which break conformal invariance "by 



hand" [39, IIO, 41, 42]. These models are also able to fit the thermodynamical functions of 
QCD at temperatures close to the phase transition. 
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Figure 12: Comparison of LO, NLO, and NNLO predictions for the scaled trace anomaly with 
SU(3) pure-glue lattice data from Boyd et al. p^. Shaded bands show the result of varying the 
renormalization scale /i by a factor of two around /i = 2t:T. 



In Fig. 13, we show the HTLpt predictions for the trace anomaly scaled by T'^/T^ using 
the BN mass prescription and three-loop running of a^. The points are lattice data from 



Boyd et al. [28|. The most remarkable feature is that lattice data are essentially constant 
over a very large temperature range. Clearly, HTLpt does not reproduce the scaled lattice 
data precisely; however, the agreement is dramatically improved when going from NLO to 
NNLO. 



24 



1.2 

Si.o 

bo.8 

«0.6 

SO.4 

Oh 
70.2 



LO HTLpt 

NLO HTLpt 

NNLO HTLpt 

Lattice Data (Boyd et al) 




0.0: 



Figure 13: Comparison of LO, NLO, and NNLO predictions for the scaled trace anomaly multiplied 
by T^/T^ with SU(3) pure-glue lattice data from Boyd et al. [|8|. Shaded bands show the result of 
varying the renormalization scale /i by a factor of two around fi = 2ttT. Orange band shown with 
lattice data indicates scaled error assuming a ±2.5% error in the lattice data for the trace anomaly. 

6. Conclusions 

In this paper, we have presented results for the LO, NLO, and NNLO thermodynamic 
functions for SU(A^c) Yang-Mills theory using HTLpt. We compared our predictions with 
lattice data for Nc = 3 and found that with a perturbative mass prescription HTLpt is 
consistent with available lattice data down to approximately T ~ 3 Tc in the case of the 
pressure and T ~ 2 Tc in the case of the energy density and entropy. These results are in line 
with expectations since below T ~ 2 — 3 Tc a simple "electric" quasiparticle approximation 
breaks down due to nonperturbative chromomagnetic effects. 

The mass parameter 771 £> in HTLpt is arbitrary and we employed two different prescrip- 
tions for fixing it. Unfortunately, the variational gap equation has four complex conjugate 
solutions, two with positive real parts. This has also been observed in scalar theory and 
QED. Whether this is a problem of HTLpt as such or is related to our m£,/T expansion 
is unknown. Since it is not currently possible to evaluate the NNLO HTLpt diagrams in 
gauge theories exactly, it is impossible to settle the issue at this stage. On the other hand, 
the BN mass prescription is well defined to all orders in perturbation theory and does a 
reasonable job reproducing available lattice data for temperatures above T ~ 3Tc. 

That being said without lattice data to compare with one would be hard pressed 
to favor one prescription over the other. While it is true that variational solutions are 
complex, one could be tempted to ignore the imaginary contributions and take Re['P] as 
the approximation. As we have shown (see Fig. |^ in this case the convergence of HTLpt 
is impressive; however, the result lies above the lattice data at low temperatures. Taking 
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o 
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o 
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o 
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the BN mass prescription on the other hand results in a real-valued pressure which seems 
to be in better agreement with the lattice data at the expense of poorer convergence of 
the successive approximations. The dependence on the mass prescription gives another 
measure of our theoretical uncertainty. If one were able to carry the HTLpt 5 expansion to 
higher and higher orders the dependence on the mass prescription would go away; however, 
going to higher orders presents a rather daunting task since one encounters a purely non- 
perturbative contribution at 0{a^) for which lattice input is required. 

Finally, we emphasize that HTLpt is gauge invariant by construction and that it can 
be used to calculate time-dependent quantities as well since it is formulated in Minkowski 
space. The fact that our BN mass results are consistent with lattice data down to approx- /- 

imately 3Tc, suggests that HTLpt can provide a coherent framework that can be used to | , 

systematically calculate real-time quantities such as heavy quark diffusion and viscosity at j~^ 

temperatures relevant for LHC. L J 
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A. HTL Feynman rules 

In this appendix, we present Feynman rules for HTL perturbation theory in QCD. We give 
explicit expressions for the propagators and for the 3- and 4-gluon vertices. The Feynman 
rules are given in Minkowski space to facilitate future application to real-time processes. A 
Minkowski momentum is denoted p = (pq, p), and the inner product is p ■ q = poQo — p • q. 
The vector that specifies the thermal rest frame is n = (1,0). 

A.l Gluon self-energy 

The HTL gluon self-energy tensor for a gluon of momentum p is 

U^'^ip) = ml [T^^ip, -p) - n^n''\ . (A.l) 

The tensor T'^'^ip, q), which is defined only for momenta that satisfy p + q = 0, is 

r'''^(p,-p) = /yV — \ . (A.2) 

\ p-y/y 

The angular brackets indicate averaging over the spatial directions of the light-like vector 
y = (l,y). The tensor T^"^ is symmetric in ^ and ly and satisfies the "Ward identity" 

Pt^T^^iP, -P) = P-n n" . (A. 3) 
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The self-energy tensor H^'^ is therefore also symmetric in /x and v and satisfies 

p^n'^^'ip) = , (A.4) 

g^.U^-'ip) = -ml . (A.5) 

The gluon self-energy tensor can be expressed in terms of two scalar functions, the 
transverse and longitudinal self-energies Ht and Tii, defined by 

^t{p) = ^{5'' -flP)li''{p) , 
liLip) = -Tf\p) , 

where p is the unit vector in the direction of p. In terms of these functions, the self-energy 
tensor is 

u'^-'ip) = -^T(p)^7-^^z.(p)L^^ 

Up 

where the tensors T„ and L„ are 






Lr 



ripnTp 



P nl 



(A.6) 




(A.7) 


Q 


energy 


m 




±j 


(A.8) 


•n 




o 




o 


(A.9) 


INJ 


(A.IO) 


O 



nf'p = n^'- -^p" , (A.ll 



{d - 1)Ut{p) + ^Ul{p) = ml . (A.12) 

np 



We can express both self-energy functions in terms of the function T^'' defined by (A. 2): 

2 

ULip)=ml[l-r'^{p,-p)], (A.14) 

In the tensor T^'^{p,—p) defined in (|A.2D , the angular brackets indicate the angular 
average over the unit vector y. In almost all previous work, the angular average in ( [A.2[ ) 
has been taken in d = 3 dimensions. For consistency of higher-order corrections, it is 
essential to take the angular average in d = 3 — 2e dimensions and analytically continue to 
d = 2> only after all poles in e have been cancelled. Expressing the angular average as an 
integral over the cosine of an angle, the expression for the 00 component of the tensor is 

T°°(P, -P) = ^fdc (1 - c2)-^^V ' (^-15) 

2 J -I Po- \p\c 
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The four-vector rip is O 



n-p 



and satisfies p-Up = and 11^ = 1 — {n-p)'^ /p^ . The equation ( A.5 ) reduces to the identity CD 

o 



where the weight function w{e) is 



He) = ££^2- 



r( 



r(i)r(i-6) 



(A.16) 



The integral in ( A.15| ) must be defined so that it is analytic at po = oo- It then has a 

IpI to po = +|p|- If we take the limit e — > 0, it reduces to 



branch cut running from po 



7^\P,-P) 



Po ■, Po+ P 
■log- 



2|p| Po - IPI 
which is the expression that appears in the usual HTL self-energy functions. 

A. 2 Gluon propagator 

The Feynman rule for the gluon propagator is 



(A.17) 



(A.18) 



where the gluon propagator tensor A^i, depends on the choice of gauge fixing. We consider 
two possibilities that introduce an arbitrary gauge parameter ^: general covariant gauge 
and general Coulomb gauge. In both cases, the inverse propagator reduces in the limit 
^ ^ oo to 



^^Hpr = -p'g'"' +p^p'' - n^"(p) . 



This can also be written 



1 1 

A'^^ivY" = — T^"" H L^"" 

"^^^P' At{p) P ^njALip) P ' 

where At and A^ are the transverse and longitudinal propagators: 

1 



At{p) 
^l{p) 



The inverse propagator for general ^ is 



p2 - IiT{p) ' 
1 



1 



^-^{pY^ = A^\pY^ - ^p^'p" covariant , 

= A^ipY" - \ [p^ -p-n n^) {p'' -p-n n^) 

Coulomb . 



The propagators obtained by inverting the tensors in ( |A.24| ) and ( [A. 23 ) are 



A^^'^ip) = -At{p)TI:'' + AL{p)n^n; - C 



T)^T) 



(p2)2 

covariant , 



-AT(p)r7 + AL(pKn''-f 



2^2 V 



[ripp 
Coulomb . 



(A.19) 
(A.20) 

(A.21) 

(A.22) 

(A.23) 
(A.24) 



(A.25) 



(A.26) 
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It is convenient to define the following combination of propagators: 



Axip) = Al(p) + ^At(p) . 



nt: 



(A.27) 



Using ( A.12| ), ( A.21| ), and (A. 22), it can be expressed in the alternative form 

Ax(p) = [ml - dUrip)] Al{p)At{p) , (A.28) 

which shows that it vanishes in the limit rriD —> 0. In the covariant gauge, the propagator 
tensor can be written 

A^-(p) = [-At(p)5^^ + Ax{p)nV] - ^Ax(p) {p^n'' + n^p") 



+ 



At(p) + ^Ax(p)-4 

pA pZ 



pfJ.pl' 



(A.29) 



This decomposition of the propagator into three terms has proved to be particularly con- 
venient for explicit calculations. For example, the first term satisfies the identity 



-AT{p)g^,u + Ax{p)nf,n^\ A^ [p) = g^ - -^ + n 2 a f \ Pf'^p 

f pz n^p'^ Ai[p) ^ 



(A.30) 



A. 3 Three-gluon vertex 



The three-gluon vertex for gluons with outgoing momenta p, q, and r, Lorentz indices //, 
u, and A, and color indices a, b, and c is 

iKbciP' (?, r) = -gfabcr''''\p, q, r) , (A.31) 

where /" are the structure constants and the three-gluon vertex tensor is 

r^""(p, g, r) = 5^-(p - qf + g''\q - vY + g^^{r - pY " mlT^''\p, q, r) . (A.32) 

The tensor T'^'^a j^^ ^.j^g HTL correction term is defined only for p + q + r = ^■. 



T>'''\p,q,r) 



y^y^y'f ^'^ 



r-n 



p-y q-y r -y q-y 



(A.33) 



This tensor is totally symmetric in its three indices and traceless in any pair of indices: 
g^j^vT^^ = 0. It is odd (even) under odd (even) permutations of the momenta p, q, and r. 
It satisfies the "Ward identity" 



Q^T^'^Hp, q, r) = T'Hp + q,r)- T^\p, r + q) . 
The three-gluon vertex tensor therefore satisfies the Ward identity 

p,T^^^\p,q,T) = A-J{qY^-A-J{rY^. 



(A.34) 



(A.35) 
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A. 4 Four-gluon vertex 

The four-gluon vertex for gluons with outgoing momenta p, q, r, and s, Lorentz indices //, 
u, A, and a, and color indices a, b, c, and d is 



+2m|,tr 



T^''^''{p,q,r,s)} 



rpa I rpbrpcrpd _i_ rpdrpCrph ^ 

+ 2 cyclic permutations , (A. 36) 

where the cyclic permutations are of {q,i^,b), (r, A,c), and {s,a,d). The matrices T°- are 

the fundamental representation of the SU{3) algebra with the standard normalization (~ i 

tr(r"r ) = 2^"'^- The tensor 7"/^'^-'*°" in the HTL correction term is defined only for p + q+ Hp i 

\ \p-y q-y {q + r)-y ^~^ 

^^P + qyn ^ {p + q + r).n \ \ Q 
q.yr-y {r + s)-y r-y s-y [s + p)-y J j 

This tensor is totally symmetric in its four indices and traceless in any pair of indices: l\^ 



r + s = 0: 



(A.39) 

When its color indices are traced in pairs, the four-gluon vertex becomes particularly 
simple: 

^"'^"'KlcdiP^ ^> ^> «) = -w''N,{Nl - l)r^-'^'^(p, q, r, s) , (A.40) 

where the color-traced four-gluon vertex tensor is 

Yi^^M{j,^ g, r, s) = 2g''^g^'' - gi'^g"'' - g'^'g"^ - m^T'"'^'' {p, s, q, r) . (A.41) 

Note the ordering of the momenta in the arguments of the tensor T'M'^-^o'^ which comes 



from the use of the Bianchi identity (A. 391). The tensor (A.41) is symmetric under the 



interchange of ^ and v, under the interchange of A and a, and under the interchange 
of (//, z^) and (A, a). It is also symmetric under the interchange of p and q, under the 
interchange of r and s, and under the interchange of {p, q) and (r, s). It satisfies the Ward 
identity 



p^r'^''^^^p,q,r,s) = T''^^q,r+p,s)-V'^'^{q,r,s+p) . (A.42) 



o 



g^ivT^^ = 0. It is even under cyclic or anti-cyclic permutations of the momenta j», q, r, 
and s. It satisfies the "Ward identity" 

q^.T'^'^'^ip, q, r, s) = T^^^ip + q, r, s) Q 

-'r^''{p,r + q,s) (A.38) ^^ 

and the "Bianchi identity" C J 

T^-^^p, q, r, s) + T^'^^'^ip, r, s, q) + T^'^^^p, s, q,r)=0 



o 
o 
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A. 5 Ghost propagator and vertex 

The ghost propagator and the ghost-gluon vertex depend on the gauge. The Feynman rule 
for the ghost propagator is 

—(5"^ covariant , (A. 43) 



P 



Coulomb . (A.44) 



The Feynman rule for the vertex in which a gluon with indices n and a interacts with an 

outgoing ghost with outgoing momentum r and color index c is (~ , 

_gfabc^f, covariant , (A. 45) PC 

-gr''''{r^-r-nn^') Coulomb. (A.46) M 

o 

A. 6 HTL counterterm C !) 

The Feynman rule for the insertion of an HTL counterterm into a gluon propagator is ^-"^ 



Every closed ghost loop requires a multiplicative factor of — 1. 



-i6''''U^"'{p) , (A.47) 

where 11^'^ (p) is the HTL gluon self-energy tensor given in (A.^ 
A. 7 Imaginary-time formalism 



p-n — > —P-N where N = (— i,0), p-y — > —P-Y where Y = (— i,y), and y^ — )• Y^. For 
example, the Euclidean tensor corresponding to ([A.2|) is 



(--^) 



T^'^iP, -P) = ( Y'^Y''—— ) . (A.48) 



The average is taken over the directions of the unit vector y. 

Alternatively, one can calculate a diagram by using the Feynman rules for Minkowski 
momenta, reducing the expressions for diagrams to scalars, and then make the appropriate 
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In the imaginary-time formalism, Minkoswski energies have discrete imaginary values po = 

i(27rnT) and integrals over Minkowski space are replaced by sum-integrals over Euclidean 

vectors (27rnT, p). We will use the notation P = {Pq,p) for Euclidean momenta. The CZD 

magnitude of the spatial momentum will be denoted p = |p|, and should not be confused ^ ^ 

with a Minkowski vector. The inner product of two Euclidean vectors is P-Q = PoQo+P'Q- 

The vector that specifies the thermal rest frame remains n = (1, 0). 

The Feynman rules for Minkowski space given above can be easily adapted to Euclidean 
space. The Euclidean tensor in a given Feynman rule is obtained from the corresponding 
Minkowski tensor with raised indices by replacing each Minkowski energy po by iPq, where 
Pq is the corresponding Euclidean energy, and multipying by — i for every index. This 
prescription transforms p = {pq, p) into P = {Pq, p), g^'^ into —6^'^, and pq into —PQ. The 
effect on the HTL tensors defined in ([A. 21) , ( [A. 331 ), and ( A.37| ) is equivalent to substituting 



substitutions, such as p^ — t- — P^, p ■ q ^- —P • Q, and n ■ p ^- in ■ P. For example, the 
propagator functions ( A.21| ) and ( [A.22 ) become 



At(P) 



-1 



P2 + Ut{P) ' 

1 



(A.49) 
(A.50) 



The expressions for the HTL self-energy functions I[t{P) and IIl{P) are given by ( A.13| ) 
and ( [A. 14 ) with ni replaced by rip = p'^/P'^ and T^^{p, —p) replaced by 



2 J_i iPo-pc 



(A.51) 



Note that this function differs by a sign from the 00 component of the Euclidean tensor 



corresponding to (|A.2D : 



T«0(p,-p) = -r°°(p,-p) 



-Tp 



(A.52) 



po-^iPo 

A more convenient form for calculating sum-integrals that involve the function Tp is 



Tp 



p2 



P^+P^C^ ,^ 



where the angular brackets represent an average over c defined by 



Jo 



(A.53) 



(A.54) 



and w{e) is given in ( |A.16 ). 



B. Sum-integrals 

In the imaginary-time formalism for thermal field theory, the 4- momentum P = {Pq,p) is 
Euclidean with P^ = Pq -l-p'^. The Euclidean energy Pq has discrete values: Pq = 2mrT for 
bosons, where n is an integer. Loop diagrams involve sums over Pq and integrals over p. 
With dimensional regularization, the integral is generalized to d = 3— 2e spatial dimensions. 
We define the dimensionally regularized sum-integral by 



i-^ 



Att 



Po=2n-KT 



S- 



2e 



P 



{2n 



i3-2e 



(B.l) 



where 3 — 2e is the dimension of space and fi is an arbitrary momentum scale. The factor 

(e'^^ /AttY is introduced so that, after minimal subtraction of the poles in e due to ultraviolet 

divergences, ^ coincides with the renormalization scale of the MS renormalization scheme. 

Below we list the sum-integrals required to complete the three loop calculation. We 



refer to Ref. [18| for details concerning the sum-integral evaluations. 



o 
o 

o 



o 

o 
o 
o 
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B.l One- loop sum-integrals 

The simple one-loop sum-integrals required in our calculations can be derived from the 
formulas 

p^"^ _/ /i \2. 2r(|+m-e)r(n-f -m + e) 
P(P2)^~V4^; r(n)r(2-2e) 

X r(l - e)C(2n - 2m - 3 + 2e)e'^-E 

rTi4+2m,— 2n/r) \l+2m— 2ra 



(B.2) 



The specific bosonic one-loop sum-integrals needed are 



vr 






p2 12 Uvrr; 



1+ (2 + 2C^)e + 0(e2^ 



^(P^ 



1 



(P2)2 (4^)2 V47rT 



M 



^p2p2 (47r)2 V47rr 



M 



2e 



2e 



C(-i) 



- + 27B + 0(e^ 



1 



+ 27s + 2 + 0(e) 



The number 71 is the first Stieltjes gamma constant defined by the equation 

C(l + z) = - + 7s - liz + 0{z^) . 



(B.3) 

(B.4) 

(B.5) 
(B.6) 

(B.7) 



We also need some more difficult one-loop sum-integrals that involve the HTL function 
defined in ( [A. 51 ). The specific bosonic sum- integrals needed are 



i^(^ri-i) 



(47r)2 V47rr/ 



TpP^ 

7^'^^^ " (4^ U^J 
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+ 27i5 + 21og2 + 0(e) 
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(B.8)^ 



2e 



21og2(i + 27sj +21og22 + ^ + 0(e) 



(B.9) 



2€ / 2 
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1 \ 4 22 

:i + 21og2)( - + 27sj -- + _log2 + 21og2 2 + (!)(e) 



(B.IO) 



B.2 Two-loop sum-integrals 

The simple two-loop sum-integrals that are needed are 



1 



PQ P^Q'R' 
1 



PQ 



PQ 



p2Q2^2 
p2Q2^4 



0{e), 

r2 / ^ N4e 1 

(47r)2 Vi^y 12 



T2 



(-^) 



ie 1 



(47r)2 \4ttTJ 6 



1 C'(-l) 

- + 10-121og2 + 4^^-^-^ + O(e) 

e C(-l) 

1 8 C 

- H h 27£; + 2—, 

e 3 ' C -1 



1) 



+ 0{e) 



(B.ll) 
(B.12) 
(B.13) 
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-PQ P^Q^r^R^ ~ (4^ \4^) 9 



1 



+ 7.521 + 0(e) 



(B.14) 
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r2 



^ \4e / 1 



^Q P2g2^4 (4^)2 



\4ttTJ 



i + !+41„g24,,+ lffci)+OM 



(B.15) 



where R = ~{P + Q) and r = |p + q|. We also need some more difficult two-loop sum- 
integrals tfiat involve the functions Tp defined in ( |A.51 ). The specific bosonic sum-integrals 
needed are 
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PQ 



PQ 



PQ 



p2Q2^: 



' ^ (47r)2 V47rr; V 48 



1 / C'i-i)\ 1 

^ + 2-121og2 + 4^ ^ ^ * 



-19.83 -h 0(e)] , 



^^ ,T, = ^(J^] 



4e 
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p2Q2^4 « (4^)2 V47rr/ V 576 
P2g2^4 /« (47r)2 UttT; V 96 
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26 24 
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92 log 2 -h 4- 



-477.7 + 0(e)] , 

1 / 8 C'(-i)A 1 

^ + ^ + 41og2 + 4-^ ^ ' ^ 
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(B.16) 
C(-l)\ 1 
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(B.17) 
(B.18) 
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B.3 Three-loop sum-integrals 

The three-loop sum-integrals needed are 

4 I = ^ rV ^ V 

rpQRP^Q^R^P + Q + R)^ 24(47r)2 \4ttTJ 



1 91 C'(-l) 



.2^141 + 0{e) 



C(-3) 
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PQRP^Q'^R^iQ-R)HR-P? 216(4^)2 Uttt; 
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1 73 12 64C'(-1) 

7^22^TT^-^^TTcFTy 



+ oie) 



(B.20) 



The three-loop sum-integrals were first calculated by Arnold and Zhai and calculational 
details can be found in Ref. |^. 

C. Three-dimensional integrals 

Dimensional regularization can be used to regularize both the ultraviolet divergences and 
infrared divergences in 3-dimensional integrals over momenta. The spatial dimension is 
generalized to d = 3 — 2e dimensions. Integrals are evaluated at a value of d for which they 
converge and then analytically continued to d = 3. We use the integration measure 
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C.l One-loop integrals 

The one-loop integral is given by 



[p^ + m 



2\n 



IE .."^Y 
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-m 
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Specifically, we need 



I'o = log(p^ + m^ 
Jp 



l + ^e + 0{e') 
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C.2 Two-loop integrals 

We also need a few two-loop integrals on the form 
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pg p^ + m"^ {q^ + m?)'^ (p + q)^ 
1 1 1 
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ipq V^ + rn? {q^ + m?) [(p + q)^j 
Specifically, we need Ji, J2, and Ki which were calculated in Refs. |^, [Til ]: 
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C.3 Three-loop integrals 
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We also need a number of three-loop integrals. The specific integrals we need are listed 
below and were calculated in Refs. |9|, |ll|. They are special cases of more general integrals 



defined in Ref. g|]. 
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- + 8 + 0(e) 



(C.12) 
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Finally, we need the combination 

1 (p - q)' 

{p^ + m?){q^ + m?){r'^ + m?) (q — r)^(r — p)^ 



pqr 

+ f il±^ I =(J^Y'^^\l + 0(e)] 

ipgr (P^ +"^^)[(r - P)^ +"T'^][(q-r)2 +m2] r2(p - q)2 \2mJ (Att)^ 

(C.24) 
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